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A general form of the fifth-order nonlinear evolution equation is consid- 
ered. Helmholtz solution of the inverse variational problem is used to derive 
conditions under which this equation admits an analytic representation. A 
Lennard type recursion operator is then employed to construct a hierarchy 
of Lagrangian equations. It is explicitly demonstrated that the constructed 
system of equations has a Lax representation and two compatible Hamilto- 
nian structures. The homogeneous balance method is used to derive analytic 
soliton solutions of the third- and fifth-order equations. 
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I. Introduction 



In recent years studies in fifth-order nonlinear evolution equations have re- 
ceived considerable attention primarily because these equations possess a host 
of connections with other integrable equations which play a role in diverse ar- 
eas of physics ranging from nonlinear optics [1] to Bose-Einstein condensation 
[2]. For example, Ozer and Doken [3] used a multiple-scale method to de- 
rive the fifth-order Korteweg-de Vries (KdV) equation from the higher-order 
nonlinear Schodinger equation. On the other hand, a similar method could 
also be used [4] to obtain the nonlinear Schodinger equation from fifth-order 
KdV flow [5], Sawada-Kotera equation [6] and Kaup-Kupershmidt equation 
[7]. 

Third-order evolution equations can often be solved either by the use of 
an inverse spectral method or by taking recourse to a simple change of vari- 
ables. This is true for both the linearly dispersive KdV equation and the 
nonlinearly dispersive Rosenau-Hymann equation [8]. In contrast, it is quite 
difficult to obtain solutions of the fifth-order equations. This might be an- 
other point of interest for recent studies [9] in these equations. 

In this work we derive the condition under which the general fifth-order 
nonlinear evolution equation 

ut = u^x + Auuzx + BuxU2x + Cu^Ux, u = x{x,t) (1) 

admits an analytic representation [10] or follows from a Lagrangian. Here A, 
B and C are constant model parameters. The subscripts on u denote par- 
tial derivatives with respect to that variable and, in particular, Unx = 
We use the fifth-order Lagrangian equation to define an integrable hierarchy. 
Further, we provide a Lax representation [5] and construct a bi-Hamiltonian 
structure [11] for the system. The Lagrangian approach to nonlinear evo- 
lution equation has two novel features. First, from the Lagrangians or La- 
grangian densities we can construct Hamiltonian densities [12] which form 
the set of involutive conserved densities of the system. Second, the expres- 
sion for the Lagrangian represents a useful basis to construct an approximate 
solution for the evolution equation [13,14]. We shall, however, use a direct 
method [15] to obtain explicit analytic soliton solutions. 

In Sec. II we deal with the inverse variational problem for (1) and derive 
relations between the model parameters for the equation to be Lagrangian. 
We then make use of an appropriate pseudo-differential operator to construct 
a hierarchy of equations and present results for the first few members of the 
hierarchy. In Sec. Ill we find their Lax representations and examine tlie bi- 
Hamiltonian structure. The results presented are expected to serve as a 
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useful test of integrability. We devote Sec. IV to present explicit solitonic 
solutions by using the homogeneous balance method (HB). We present some 
concluding remarks in Sec.V. 

II. Lagrangian system of equations 

In the calculus of variation one is concerned with two types of problems, 
namely, the direct and the inverse problem of Newtonian mechanics. The 
direct problem is essentially the conventional one in which one first assigns a 
Lagrangian and then computes the equations of motion through Lagrange's 
equations. As opposed to this, the inverse problems begins with the equation 
of motion and then constructs a Lagrangian consistent with the variational 
principle [10]. The inverse problem of the calculus of variation was solved by 
Helmholtz [16] during the end of the nineteenth century. For continuum me- 
chanics the Helmholtz version of the inverse problem proceeds by considering 
an r-tuple of differentiable functions written as 

P[v]^ P{x,v^''^) e A'' (2) 

and then defining the so-called Frechet derivative. The Frechet derivative of 
P is the differential operator Dp : A'^ ^ and is given by 

^P[v+eQ[v]] (3) 

for any Q e A^. The Helmholtz condition asserts that P is the Euler-Lagrange 
expression for some variational problem iff Dp is self-adjoint. When self- 
adjointness is guaranteed, a Lagrangian density for P can be explicitly con- 
structed using the homotopy formula 

C[v] = [ vP[Xv]d\ . (4) 
Jo 

In the following we shall demand the Helmholtz condition to be valid for (1). 
This will provide us with certain constraints between the model parameters 
of (1) to follow from a Lagrangian density. 

A single evolution equation Ut = P[u], -u e IR is never the Euler-Lagrange 
equation of a variational problem [16]. One common trick to put a single 
evolution equation into a variational form is to replace u by a potential 
function 

u — —Wx, w — w{x, t) . (5) 

The function w is often called the Casimir potential. In terms of the Casimir 
potential, (1) reads 

Wxt = P[wx] , (6) 



Dp{Q) = 



de 
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where 

P[w^] = Wex - Aw^Wi^ - Bw2xW3x + CwJ^W2x ■ (7) 

From (3) and (7) we obtain 

Dp = Dex - AwxD^x - Aw^xD^ - Bw2xD3x - Bw3xD2x 

+CWx'^D2x + 2CWxW2xDx . (8) 

To construct the adjoint operator D* of the above Frechet derivative we 
rewrite (8) as 

j 

and make use of the definition [16] 

D*p^J2(-^)i-Pi (10) 

meaning that for any Q e A 

j 

This gives 

D*p = Dqx - AwxD^x - {SA - B)w4xDx - (4A - B)w2xD3x 

-(6A - 2B)w3xD2x + CwlD2x + 2CwxW2xDx . (12) 

Demanding variational self-adjointness we obtain from (8) and (12) 

B^2A (13) 

while C remains unrestricted. Thus the nonlinear equation 

Ut = + AuUsx + 2AUxU2x + Cu^Ux (14) 

forms a Lagrangian system. We note that the Lax equation [5] with A = 10, 
S = 20 and C = 30 and the Ito equation [17] with A = 3, B = 6 and C = 2 
are of the form (14) while the Sawada-Kotera equation with A = B = C = 5 
and the Kaup-Kupershmidt equation with A — 10, B — 25 and C = 20 are 
non-Lagrangian. 

We now use the fifth-order Lagrangian equation (14) to define an inte- 
grable hierarchy. To that end we introduce a pseudo-differential or integro- 
differential operator A which acts on a generic function f{x) to give [18] 

r+oo 

A/(x) = fxx - puf{x) + qux / dyf{y). (15) 

J X 
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Further, we introduce a function g^'^ to follow from 

N^u,{x,t) ^gt\n ^Q,l,2 (16) 

Here is a polynomial in u and its x-derivatives (up to derivative of order 
2n). Using /(x) = Ux{x,t) in (15) we have 

A/(x) = (U2. - ^-^u% . (17) 

From (16) and (17) 

A^Ux{x,t) =U5x-{2p + q)uu3x-{3p + 4:q)uxU2x + {p + q)ip+-)u^Ux ■ (18) 

Comparing (14) and (18) and identifying A^Ux{x,t) as Uf we can express p 
and q in terms of A. This allows us to write 

C = {p + q){p + l) = -^ . (19) 

Therefore, the general form of the fifth-order Lagrangian equation generated 
by A via (16) has the form 

Ut = + AuUsx + 2AUxU2x + -^u^Ux . (20) 

We have used (16) to generate a hierarchy of nonlinear evolution equations 
for n = 0, 1,2,3 etc. The first member of the hierarchy {n = 0) is a linear 
equation given by 

Ut^Ux , (21) 
while the second one (n—1) is a third-order nonlinear equation 

2>A 

Ut = u^x + -r'^'^x ■ (22) 
o 

The third member (n = 2) is obviously the fifth-order equation given in (20) . 
The corresponding seventh and ninth order equations are given by 

7 A 21A 7^2 2 

Ut = U7x + -rUbxU + —^U^xUx + "Au3xU2x + -TTtU U^x 

5 5 10 

14^' 3 7yl3 „ 

H — -^uUxU2x + -^"^x + -^u Ux (23) 
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and 



9A 36A 84A 126A 651A'^ ^ 

Ut = U9x + -^U-jxU + —^UqxUx + —^U5xU2x H ^U^xUsx H ^^^x^a^ 

483^2 2 63^2 378^2 53^2 
H — g^-Ux^sx H ^uu2xU3x H g^^i^tx^^4x + ^^"^ 

63^3 . 126^3 2 21A'^ o 63^4 . , , 

III. Lax representation and bi-Hamiltonian structure 

Integrable nonlinear evolution equations admit zero curvature or Lax rep- 
resentation [5]. These equations are characterized by an infinite number of 
conserved densities which are in involution. Moreover, each number of the hi- 
erarchy has a bi-Hamiltonian structure [11]. In the following we demonstrate 
these three important features for our equations in (20)- (24). 

The Lax representation of nonlinear evolution equations is based on the 
algebra of differential operators. Here one considers two linear operators L 
and M . The eigenvalue equation for the operator L is given by 

Lijj = Xip (25) 

with i/j, the eigenfunction and A, the corresponding eigenvalue. The operator 
M characterizes the change of eigenfunctions with the parameter t which, in 
a nonlinear evolution equation, usually corresponds to the time . The general 
form of this equation is 

i/jt^ Mip . (26) 

If we now invoke a basic result of the inverse spectral method that ^ = 
for non-zero eigenfunctions [19], then (25) and (26) will immediately give 

^-[M,L]. (27) 

Equation (27) is called the Lax equation and L and M are called the Lax 
pairs. In the context of Lax's method it is often said that L defines the orig- 
inal spectral problem while M represents an auxiliary spectral problem. For 
a given nonlinear evolution equation one needs to find these operators. This 
is not always a straightforward task. In fact, no systematic procedure has 
been derived to determine whether a nonlinear partial differential equation 
can be represented in the form (27). 

We shall now find the Lax representation for the hierarchy of equations 
given in (20)-(24). We first note that as one goes along the hierarchy the orig- 
inal spectral problem remains invariant while the auxiliary spectral problem 
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goes on changing. Keeping this in mind we take 



L^dl + ^u (28) 

In writing (28) we have exploited the similarity between (22) and the KdV 
equation. As regards the auxiliary spectral problem we postulate that for an 
evolution equation of the form Ut = K[u] the terms in the Frechet derivative 
of K[u\ contribute additively with unequal weights to form the operator M 
such that L and M via (22) reproduces K[v] . Of course, there should not be 
any inconsistency in determining the values of the weight factors. For (22) 
the Frechet derivative of K[u] can be obtained as 

3A 

Dp = dl + —{ud, + u,). (29) 
We shall, therefore, write 

3A 

Ms = + —{bud^ + cu^) . (30) 
5 

Here the subscript 3 on M indicates that (30) represents the second Lax oper- 
ator for the third-order equation. We shall follow this convention throughout. 
Equations (22), (27), (28) and (30) can be combined to get a = 4, 6 = 1 and 
c = |. Thus we have 

M^^4dl + —{ud, + -u,) . (31) 
o z 

Similarly, we find the results 

Ms = 16dl + 4Audl + GAu^dl + bAu2A + + ^^*3x + -[^^^x ,(32) 

112 A IAA^ 
My = 649^ + —^udl + b&Au^dl + MAu2:cdl + -^u^dl + IQAu^^dl 

^"^■^ ^2 IQIA ^ ^ 147^2 2^ 7A3 3^ 
-Mw^d^ H ^u^xOx + (A uu2xOx H 3o~^^^^ + 

21A^ 21A^ 21A^ 2 / N 

■U5x + —rp-UU3x H ^UxU2x + -TT^W Ux , (33) 



42^2 
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576A ^7 51A 2016A . 897A ^ 4368A ^5 
Mg = 2569^ + —^udl + -^u^x + —^Uxdl + —^uqxOx + —^U2xdl 

2814A , 5124A ^5 5061^2 

+ «5a.5' + 1176AM3c.a^ + —^UAxdl + 252A\^M2a.5' + —^Ul^dx 



iU 
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H ^ — u^o^ H — u^usA H ^ — uu2xd^ + 126A uus^d^ 

504^2 557^2 5^g^2 2967^2 



609^2 504^2 557^2 5^g^2 2967^^ 

«2x^3x H ^^^V, d;, H ^^tiW5x H ^ MWxC'a; ^^^4x4 

1743^2 189^3 44;^^^ 3^ 42^^ 21^3 „ 

H ul H uu^Ox H H u u^x 

50 100 50 25 100 

r8A3 2 ^2 63^3 2 ^ 189 63^4 63 o 

U UxO^ H M UoxOx H UUxUox H 'i^ C'a; H u u .(34) 

50 10 25 1000 500 ^ ^ 



Zakharov and Faddeev [20] developed the Hamiltonian approach to in- 
tegrability of nonlinear evolution equations in one spatial and one temporal 
(1+1) dimension and, in particular, Gardner [21] interpreted the KdV equa- 
tion as a completely integrable Hamiltonian system with d.j. as the relevant 
Hamiltonian operator. A significant development in the Hamiltonian the- 
ory is due to Magri [11] who realized that integrable Hamiltonian systems 
have an additional structure. They are bi-Hamiltonian i.e. they are Hamil- 
tonian with respect to two diff'erent compatible Hamiltonian operators. The 
bi-Hamiltonian structure of the integrable equation is based on a mathemat- 
ical formulation that does not make explicit reference to the Lagrangian of 
the equations in the hierarchy [22]. Here we shall demonstrate that the bi- 
Hamiltonian structure of the system of equations (20)- (24) can be realized in 
terms of a set of Hamiltonian densities obtained from the Lagrangians. Using 
(4) we can obtain the Lagrangian densities for our equations. In particular, 
we have ^ ^ 

>Ci = -WtWx - -wl , (35) 

1 1 A o 

A = -WtWx - -WxWsx + —w^ , (36) 

^5 = ^WtWx - ^WxW5x + ^wlw3x + ^WxWj^ - , (37) 



1 . 7A , 7^2 , , 7^2 

40 

7A' 



and 



1 1 3y4 2 o 2 7^2 , 

£9 = -i^WtWx - -WxWqx - -yW^^W^x + -y'^Sx - -^^'^xW^^ + ^«^2x 

-Wx^sx - -T7^Wx'^2x - TF^^x- (39) 



200 100 10000 
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In the above L\ is the Lagrangian density for the linear equation in (21). The 
other subscripts on L are self explanatory. The corresponding Hamiltonian 
densities are given by 

Til = , (40) 
1 A 

= l^UU2^ + —u^ , (41) 

1 A A A?' 

'H5 = -^^4^ + —u^U2x + -^uul + — -u^ , (42) 

= ^MM6a; + ^mmL + "^^^^^ + "^^^^2x + > (43) 



and 



Tig = 2^^8x - -Y'^^'^^x + -^u^^ - - ^2x - -^^^ 

+ 10000" ■ 

In the theory of Zakharov and Faddeev [20] and of Gardner [21] the 
Hamiltonian form of an integrable nonlinear evolution equation reads 

ut = a,(^) (45) 

with 7i, the Hamiltonian densities of that equation. Here ^ denotes the 
usual variational derivative written as 

i; = Ew^. ». = ar». (46) 

n>0 " 

Using the Hamiltonian densities in (40)- (44), one can easily verify the 
Faddeev-Zakharov-Gardner equation in (45) to yield the appropriate nonlin- 
ear equations in (20)- (24). The bi-Hamiltonian form of evolution equations 
is given by [11] 

= ^^^^^^ = ^(^^ 

with m = 2n+l, n = 0, 1, 2, In (47) the second Hamiltonian operator 

is related to the recursion operator by [16] 

£ = (48) 

From (15) and (48) we get 

S = dl + ^ud^ + 4«a. • (49) 
o o 
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From (47) and (49) we have 



Wt = 9x(^^) = (<9^ + ^M9a; + -M^)(-^) . (oO) 



For n—1 (50) reads 

= 5.(^) = (9| + '^-^ud^ + y^^.)(^) . (51) 

From (41), (42) and (51) one can easily obtain (20) verifying the bi- 
Hamiltonian structure. Similar results can also be checked for other pairs 
of the Hamiltonians in (40)-(44). 

IV. Soliton solution 

We have just seen that the bi-Hamiltonian form (51) corresponds to the 
fifth-order nonlinear equation in (20). Here we shall make use of the homoge- 
neous balance method (HB) [15] to construct an analytical expression for the 
soliton solution of this equation. According to HB method, the field variable 
is first expanded as 

N 

where the superscript {%) denotes the derivative index. In particular, /"^^^ = 
If) Z'^^'' = and so on. Substituting (52) in (20) and balancing the con- 
tribution of the linear term with that of the nonlinear terms, the expression 
in (52) becomes restricted to 

i.(x,i) = /(^)^^ + /««;2. , (53) 

where subscripts on w stand for appropriate partial derivative. From (53) 
and (20) we have 



<x,t) = ^/»(^a;,t)) , (52) 



■pother terms involving lower powers of the partial derivatives of w = . 
Setting the coefficient of wl to zero we get 

fi7) + ^J(2) J(5) + 2^^(3)^(4) ^ ^(/(2))2;(3) ^ q . (55) 

If we try a solution of (55) in the form 

f — a Inw (56) 
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we immediately get 

20 , , 

a = -. (57) 

From (56) we can deduce the following results 

m f(5) ^ _^ f{i) fm m ^ _^ nr) m\2 m ^ -.(7). 

J J 20 ' 30 1 J J 5 

C /■(2)N3 ^ ^(6) m f{2) r{3) ^ ^ .(6). 

\j ) 120 ' ^ ^ 60 ' 

f (2) WS) ^ /■(5) /.(I) /■(4) ^ _« .(5) / ^(2)\2 .(1) ^ tt^ .(5) 

J J 12 1 J J 1 \j ) J 24 ' 

/(3) ^ _^ J(4)^ (J(2)^2 ^ _« y(4)^ ^y(l)^2 y(2) ^ ^ j(4). 

3 6 6 

(/(1))2 = _«y(2) _ (58) 

Substituting (58) in the full form of (54), the latter is reduced to a linear 
polynomial in /*^^\ /*^^\ /^^^ . If the coefficient of each is set equal to 
zero we get a set of partial differential equations for w{x, t) 

Wxxt - wtx = , (59a) 

2wxWxt + WtWxt+ {2Aa-?>^)wzxWAx + {Aa-2l)w2xW^x-^WxW&x = , (596) 

2wtw'l + {Aa - 4:2)wlw^x + {HAa - 2VS)WxW2xWax + {^Aa - l4Q)wxwl^ 
3^2 

+(16Aa - -^a^ - 210)w^^W3c, = , (59c) 

9^2 3y^2 

(48Aa - 630)u;^w^^ + (78Aa - 12m)wlw2xW3x 

10 5 

+(9Aq; - 210)«;^W4x = , (59rf) 

12^4^ 3^2 

(174yla —a^ - 2h2Q)wlwl^ + (48Aa - -— - 840)w>3^ = , (59e) 

5 10 

3^2 

(24y4a - -^a^ - 360)«;>2a. = (59/) 

11 



and 

{24Aa - -j^a^ - 360)wl = . {59g) 

Equation (59a) is a linear partial differential equation and can be con- 
verted to an ordinary differential equation by substituting 

w{x, t) — g{x + vt) — g{z) . (60) 

Using (60) in (59a) we have 

Here v is the velocity of the travelling wave represented by w{x,t). Equation 
(61) can be solved to write 

w{x, t) = g{x + vt) = Co + cie^(^+''*) , (62) 

where Cq and Ci are arbitrary constants. Using (56), (57) and (62) in (53) we 
get the exact soliton solution of the fifth-order equation in (20) and/or (51) 
in the form 

A similar result for the third-order equation in (22) is given by 

20 coCit;e^("+^*) 
"^("'^^^ A(co + c,eV^(^+-))^- ^^'^ 

The subscripts on u{x, t) are self explanatory. It is of interest to note that 
for ci = Co = 1, ^ = 10 and v = 4ac^, U3{x.t) in (64) becomes 

U3(x, t) = 2K^sech^{KX + AkH) . (65) 

From the inverse spectral method [23] for solving the KdV equation, we know 
that has a simple physical meaning. For example —n"^ represents a discrete 
energy eigenvalue of the Schodinger equation for the initial potential ^3(0;, 0). 
As in ref. 9 we shall now examine the spatial behaviour of ?i5(.x, f) at t = 0. 
For the sake of simplicity we shall work with v = 1. In Fig.l we plot ^5(0;, 0) 
as function of x for different values of parameter co and ci. All the curves 
in the figure are of sech? shape indicating that solution obtained from (63) 
have indeed solitary wave properties. The solid curve for Co = 1 and Ci = 1 
is centred at the point x = 0. If cq and ci are made unequal, the centre of 
the soliton moves either to the left or to the right. In particular, for Co > Ci, 
the shift of the centre is towards the right and we have a reverse situation 
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X 



Figure 1: Variation of U5{x,0) with x 

for Co < Ci. We have displayed this property by using a dashed curve with 
cross (cq = 4 and Ci = 1) and a simple dashed curve (cq = 1 and Ci = 4). 

V. Conclusion 

Fifth-order nonlinear evolution equations, on the one hand, have a host 
of connections with other important integrable equations and, on the other 
hand, can not be solved by simple analytical methods. These two points in- 
spired us to construct a general fifth-order equation which follow from a La- 
grangian. It is often desirable that equations of mathematical physics should 
be derivable from an action principle because a non-Lagrangian system does 
not allow one to carry out a linear stability check [24] as well as to derive a 
field theory [25] for particles described by its solutions. The Lagrangian ap- 
proach to nonlinear evolution equations is quite interesting because here one 
can derive all physico-mathematical results from first principles [8]. Based 
on the fifth-order Lagrangian equation we derived an integrable hierarchy. 
As a test of integrability we provided a Lax representation and constructed 
two compatible Hamiltoinan structures. 

We treated the third- and fifth-order equations in the hierarchy by the 
homogeneous balance method [15] to obtain analytical results for soliton 
solutions. Ideally, we could have tried the bi-linear method of Hirota [26] 
to deal with the problem because this method is very convenient for find- 
ing single- and multi-soliton solutions of nonlinear evolution equations. For 
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higher-order equations, Hirota transformation often leads to multihnear rep- 
resentation [27]. This tends to pose problems in solving the equations. The 
homogeneous balance method, on the other hand, does not involve any new 
mathematical complication as one moves from lower- to higher-order equa- 
tions. Admittedly, the algebra become more and more involved as we go up 
the ladder inside the hierarchy. The symbolic computations facilities like the 
Mapple and Mathematica can be used to circumvent algebraic complications. 



Acknowledgement. This work is supported by the University Grants Com- 
mission, Government of India, through grant No. F.10-10/2003(SR). 



14 



[I] K. Nakkeeran, K. Porsezian, P. Shanmugha Sundaram and A. Ma- 
halingam, Phys. Rev. Lett. 80, 1425 (1998). 

[2] U. Al Khawaja, H. T. C. Stoof, R. G. Hulet, K. E. Strecker, and G. 
B. Patridge, Phys. Rev. Lett. 89, 200404 (2002). 

[3] M. N. Ozer and F. T. Doken, J. Phys. A.: Math. Gen. 36, 2319 
(2003). 

[4] M. N. Ozer and I. Dag, Hadronic. J. 24, 195 (2001). 

[5] P. D. Lax, Comm. Pure AppL Math. 21, 467 (1986). 

[6] K. Sawada and T. Kotera, Prog. Theo. Phys. 51, 1355 (1974). 

[7] D. J. Kaup, Stud.AppkMath. 62, 189 (1980); B. Kupersmidt, Com- 
mun. Math. Phys. 99, 51 (1988). 

[8] P. Rosenau and J. M. Hymann, Phys. Rev. Lett. 70, 564 (1993); B. 
Talukdar, J. Shamanna and S. Ghosh, Pramana. J. Phys. 61, 99 (2003); S. 
Ghosh, U. Das and B. Talukdar, Int. J. Theor. Phys. 44, 363 (2005) 

[9] W. Hong and Y. Jang, Z. Naturforsch. 54a, 549 (1999). 

[10] R. M. SantiU, Foundations of Theoretical Mechanics I (Springer- 
Verlag, New York, 1984). 

[II] F. Magri, J. Math. Phys. 19, 1548 (1971). 

[12] B. Talukdar, S. Ghosh, J. Shamanna and P. Sarkar, Eur. Phys. J. 
D21, 105 (2002); B. Talukdar. S. Ghosh and U. Das, J. Math. Phys. 46, 
043506 (2005). 

[13] D. Anderson, Phys. Rev. A. 27, 3135 (1983). 

[14] F. Cooper, C. Lucheroni, H. Shepard and P. Sodano, 
[ar'xiv:hep-ph/9210226. vl, 9 Oct 92. 

[15] M. Wang, Phys. Lett. A. 199, 169 (1995); ibid 213, 279 (1996); 
216, 67 (1996). 



15 



[16] P. J. Olver, Application of Lie Groups to Differential Equation, 
(Springer- Verlag, NY, 1993). 

[17] M. Ito, J. Phys. Soc. Japan 49, 771 (1980). 

[18] F. Calogero and A. Degasperis, Spectral Transform and Soliton 
(North-Holland Publising Company, New York, 1982). 

[19] K. Chadan and P. C. Sabatier, Inverse problems in Quantum Scat- 
tering Theory (2nd ed. Springer, New York, 1989). 

[20] V. E. Zakharov and L. D. Faddeev, Funct. Anal. Phys. 5, 18 (1971). 

[21] C. S. Gardner, J. Math. Phys. 12, 1548 (1971). 

[22] S. Ghosh, B. Talukdar and J. Shamanna, Czech. J. Phys. 53, 425 
(2003). 

[23] C. S. Gardner, J. M. Greene, M. D. Kruskal and R. M. Miura, Phys. 
Rev. Lett. 19, 1095 (1967). 

[24] B. Dey and A. Khare, J. Phys. 33, 5335 (2000). 

[25] S. A. Hojman and L. C. Shepley, J. Math. Phys. 32, 142 (1990). 

[26] R. Hirota, Phys. Rev. Lett. 27, 1192 (1971). 

[27] J. Hietarinta in Nonlinear Dynamics: Integrability and Chaos ed. M. 
Daniel, K. M. Tamizhmani and R. Sahadevan (Narosa Publising House, New 
Delhi (2000)). 



16 



